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ABSTRACT 

It has recently been proposed that one sub-class of type la supernovae (SNe la) is sufficiently both 
distinct and common to be classified separately from the bulk of SNe la, with a suggested class name 
of "type lax supernovae" (SNe lax) , after SN 2002cx. We show that the population properties of this 
class can be understood if the events originate from helium double-detonation sub-Chandrasckhar 
mass explosions, in which a carbon-oxygen white dwarf (CO WD) accumulates a helium layer from 
a non-degenerate helium star. We have incorporated detailed binary evolution calculations for the 
progenitor systems into a binary population synthesis model to obtain birthrates and delay times for 
such events. The predicted Galactic event rate is ~0.6 — 1.8 x 10~ 3 yr" 1 , in good agreement with the 
measured rates of SNe lax. In addition, predicted delay times are ~70Myr— 710 Myr, consistent with 
the fact that SNe lax have so far only been discovered in late-type galaxies. Based on the CO WD 
mass at explosion and previous detonation models, we also estimate the distribution of resulting SN 
brightness (—13 > Mbo\ ^ —19 mag), which can reproduce the empirical diversity of SNe lax. We 
speculate on why binaries with non-degenerate donor stars might lead to SNe lax if similar systems 
with degenerate donors do not and suggest that the higher mass of the helium layer necessary for 
ignition at the lower accretion rates typically delivered from non-degenerate donors might provide the 
explanation. 

Subject headings: binaries: close — stars: evolution — supernovae: general 



1. INTRODUCTION 

Type lax supernovae (SNe lax) have been proposed to 
form a distinct sub-class of sub-luminous SNe la, contain- 
ing SNe resembling the prototype event SN 2002cx (Li et 
al. 2003; Foley et al. 2012). Those SNe lax are spectro- 
scopically similar to SNe la, but have lower maximum- 
light velocities (2000 < \v\ < 8000 km s" 1 ), typically 
lower peak magnitudes (—14.2 > My > —18.9 mag)p 
and maximum-light spectra that typically resemble those 
of the bright 1991T-likc events. Since the estimated rate 
of SN lax is roughly one third of the SN la rate, they are 
relatively common astrophysical events, although only 25 
members of the class are currently identified (Foley et al. 
2012). 

SNe lax appear not to obey the standard luminosity- 
width relation of SNe la (Foley et al. 2012 and refer- 
ences therein) ; clearly this would affect any use of them 
as distance indicators. However, perhaps more impor- 
tantly, this behaviour of SNe lax also gives us an opportu- 
nity to help us understand the physics of thermonuclear 
supernovae in general, since whatever mechanism pro- 
duces SNe lax leads to an alternative family of lightcurve 
shapes from standard SNe la. Deducing the progenitors 
of these explosions should help us to understand how 
this family of explosions differs from standard SNe la. 
Observations support the supposition that SNe lax are 
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from thermonuclear explosions of carbon-oxygen white 
dwarfs (CO WDs), due to the evidence of C/O burning in 
their maximum-light spectra (Foley et al. 2012). SNe la 
lack helium in their spectra, yet two SNe lax do show 
strong helium lines in their spectra, and so there must 
be helium in their progenitor systems. However, there 
is no hydrogen in any SN lax spectra, and significantly 
less hydrogen than helium is typically required to cause 
a signature in an SN spectra (see, e.g., Hachinger et al. 
2012). That spectral evidence suggests that in the SN lax 
progenitor systems a CO WD is accreting from a non- 
degenerate helium star or a He WD. Of those, the CO 
WD + He WD systems would exist in old stellar popula- 
tions, which is inconsistent with the observation that so 
far SNe lax have exclusively been discovered in late-type 
galaxies (e.g., Valenti et al. 2009; Foley et al. 2012). How- 
ever, CO WD + He star systems are expected primarily 
in young stellar populations, as has been observed. In 
this Letter, we investigate the population properties of 
CO WD + He star systems and ask whether they are 
consistent with being the progenitors of SNe lax. 

A CO WD can accrete material from a helium star 
to increase its mass until it ignites near to the Chan- 
drasekhar mass limit (e.g., Wang et al. 2009a). However, 
standard Chandrasekhar mass explosion models have dif- 
ficulty in reproducing the low luminosities of SNe lax 
(Hillebrandt & Niemeyer 2000). It is therefore natu- 
ral to consider whether SNe lax might be produced by 
sub-Chandrasckhar mass explosions, in which the explo- 
sion of a CO WD is triggered by the detonation of a 
substantial surface layer of accreted helium (see, e.g., 
Nomoto 1982; Woosley et al. 1986). However, the details 
of what happens following burning of the helium layer 
are still unclear. If helium ignites at the bottom of the 
helium layer, this may result in an event known as helium 
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double-detonation, during which one detonation wave 
propagates outward through the helium layer, whereas 
an inward propagating pressure wave compresses the CO 
core and leads to ignition followed by an outward deto- 
nation (see, e.g., Livne 1990; Woosley & Weaver 1994). 
Helium double-detonation sub-Chandrasekhar mass ex- 
plosions have previously been considered as promising 
explanations for standard SNe la (see, e.g., Branch et al. 
1995; Hoflich & Khokhlov 1996), but modern models dis- 
agree over which known SNe, if any, such systems might 
produce (see, e.g., Fink et al. 2007, 2010; Kromer et al. 
2010; Woosley & Kasen 2011). 

There are numerous complications, including that the 
properties of the helium layer that lead to helium igni- 
tion itself are expected to be a function of CO WD mass 
and accretion rate (see, e.g., Bildsten et al. 2007; Shen 
& Bildsten 2009). Heating from differential rotation in 
the accreted layer provides further uncertainty (Yoon & 
Langer 2004). Transients have also been observed which 
are consistent with the detonation of a thick helium layer 
on a WD which does not lead to a double-detonation 
(Poznanski ct al. 2010; Kasliwal et al. 2010). Previ- 
ous studies have often made the simplification that, for 
sufficiently massive CO WDs, a helium layer with mass 
0.1 Mq can ignite and lead to a double detonation (e.g., 
Ivanova & Taam 2004; Ruiter et al. 2011), and Fink et 
al. (2007, 2010) support the position that layers of that 
mass can produce a double-detonation. 

Generic double-detonation scenarios have previously 
been suggested for SNe lax, but inconclusively (see, e.g., 
Foley et al. 2012 and references therein). We note that 
Foley et al. (2009) suggested that an alternative explo- 
sion model might be responsible for SNe lax, specifically 
a failed deflagration model (see, e.g., Foley et al. 2012). 
In models like this, the accreting WD would survive and 
potentially possess peculiar observational properties (see 
also Jordan ct al. 2012; Kromer et al. 2012). 

Overall, determining the progenitors of SNe lax could 
distinguish between theoretical models which are largely 
separated by precise details of thermonuclear explosion 
physics at WD densities. 

In U2]we further describe our assumptions and binary 
evolution calculations. Our population synthesis results 
are presented in §3\ followed by a discussion and sum- 
mary in §4. 

2. BINARY EVOLUTION CALCULATIONS 

Employing Egglcton's stellar evolution code (Egglcton 
1973; later updated by Han et al. 1994; Pols et al. 1998), 
we have calculated the evolution of the CO WD + He 
star systems. Roche lobe overflow (RLOF) is treated 
within the code as described by Han et al. (2000). We set 
the ratio of mixing length to local pressure scale height, 
a = I /Hp, to be 2.0. In our calculations, the initial 
helium star models are composed of helium abundance 
Y = 0.98 and metallicity Z = 0.02. Orbital angular mo- 
mentum loss due to gravitational wave radiation (GWR) 
is included. 

In the double-detonation model, the helium star trans- 
fers some of its material onto the surface of the WD, 
which increases the mass of the WD as a consequence. 
If the mass-transfer rate onto the WD from the helium 
star is higher than 4 x 10 -8 Mq yr _1 , stable burning al- 
lows the CO mass of the WD to increase (Woosley et al. 



1986; see also Wang et al. 2009a and references therein). 
For lower accretion rates (1 x 10 _9 A/ Q yr _1 < \J\d2\ < 
4 x 10~ 8 Afoyr^ 1 ), a thick layer of helium is believed to 
grow on the surface of the WD. When the mass-transfer 
rate drops even further (\M 2 \ < 1 x 10~ 9 A/ Q yr _1 ), the 
flash when the helium layer ignites has been suggested 
to be too weak to initiate a carbon detonation, which re- 
sults in only a single helium detonation wave propagating 
outward (see, e.g., Nomoto et al. 1982)0 For accretion 
rates of a few times 10 -8 M Q yr _1 , Yoon & Langer (2004) 
have argued that heating by frictional dissipation signif- 
icantly reduces the eventual chance of a helium detona- 
tion, leading to some uncertainty in these mass-transfer 
rate boundaries. 

According to recent hydrodynamic simulations, the 
minimum WD mass for carbon burning might be 
~0.8M Q , since the detonation of the CO WD may be 
not triggered for lower mass (e.g., Sim et al. 2012). We 
also expect that the initial CO WD masses are below 
~1.1 Mq as more massive WDs - at formation - usually 
consist of oxygen and neon (i.e., ONe WDs). In prin- 
ciple, the WD could increase its CO mass by accreting 
helium as long as the mass-transfer rate is higher than 
4 x 10~ 8 Mq yr _1 , but this rarely occurs in our calcula- 
tions (only helium donors with masses above Mq 
lead to an increase in CO mass). Following the previous 
work described in the introduction, we assume that a 
double-detonation occurs when a helium layer with mass 
0.1 Mq accumulates on the surface of the WD. 

We incorporated the prescriptions above into Egglc- 
ton's stellar evolution code and followed the evolution of 
an ensemble of CO WD + He star systems. The mass 
lost from these systems is assumed to take away the spe- 
cific orbital angular momentum of the accreting WD. We 
have calculated the evolution of about 600 WD + He star 
systems, thereby obtaining a large, dense model grid. 
The initial mass of the helium donor stars, M%, ranges 
from 0.3 M© to 1.3 M©; the initial mass of the CO WDs, 
MJ^p, is from 0.8 Mq to 1.10 Mq; the initial orbital pe- 
riod of the binary systems, P 1 , changes from the mini- 
mum value, at which a helium zero-age main-sequence 
(ZAMS) star would fill its Roche lobe, to -0.2 d, where 
the helium star fills its Roche lobe at the end of the he- 
lium MS. 

In Figure 1, we present an example of binary evolution 
calculations for the double-detonation model. The left 
panel shows the mass-transfer rate and the mass of the 
WD envelope varying with time after the helium star fills 
its Roche lobe, whereas the right panel is the evolution- 
ary track of the helium donor star in the Hertzsprung- 
Russell diagram, where the evolution of the orbital pe- 
riod is also shown. The WD + He star binary starts with 
{Ml M^ D , log(Pyday)) = (0.60, 1.0, -1.4), where M|, 
M(y D are the initial masses of the helium star and of 
the CO WD in solar mass, and the P 1 is the initial or- 
bital period in days. Due to the short initial orbital pe- 
riod (0.04 d), angular momentum loss induced by GWR 

5 Our upper accretion rate limit only directly affects a small 
fraction of our population. The lower limit is uncertain but might 
approximately be justified in an additional way, i.e., that below 
such rates the helium layer mass needed for ignition is too large to 
normally be reached (see, e.g., Shen & Bildsten 2009). 
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Figure 1. A representative example of binary evolution calculations. Left panel: the solid and dash-dotted curves show the mass-transfer 
rate and the mass of the helium layer on the WD varying with time after the helium star fills its Roche lobe, respectively. Right panel: the 
evolutionary track of the donor star is shown as a solid curve and the evolution of orbital period is shown as a dash-dotted curve. Dotted 
vertical lines in both panels indicate the position where the doublc-dctonation may happen. 

whereas systems beyond the right boundary experience 
mass-transfer when the helium star evolves to the sub- 
giant stage. The upper boundaries are set mainly by a 
high mass-transfer rate due to orbit decay induced by 
GWR and large mass-ratio, which makes the WD grow 
in mass to the Chandrasekhar mass limit (see Wang et 
al. 2009a). The lower boundaries are where the mass- 
transfer rate M2 is higher than 1 x 10 -9 AfQyr -1 for just 
long enough to produce a 0.1 M© helium layer on the CO 
WD. 

3. BINARY POPULATION SYNTHESIS 

In this double-detonation model, the progenitor sys- 
tems containing a CO WD + He star in a close binary 
have most likely emerged from the common-envelope 
(CE) evolution of a giant binary system. CE ejection is 
still an open problem. We use the standard energy equa- 
tions to calculate the output of the CE phase (following 
Webbink 1984), in which there are two uncertain param- 
eters, a CG (the CE energetic ejection efficiency) and A 
(a structure parameter that depends on the evolutionary 
stage of the donor star and the definition of the core- 
envelope boundary). As in previous studies (e.g., Wang 
et al. 2009b), we combine a ce and A into a single free pa- 
rameter a C0 A, and show results for two values: 0.5 and 
1.5. 

In order to obtain event rates and delay times for the 
double-detonation model, we performed a series of Monte 
Carlo BPS simulations. As initial conditions we assumed 
the initial mass function of Miller & Scalo (1979) and a 
constant initial mass-ratio distribution. The distribution 
of separations is also taken to be constant in log a for 
wide binaries, where a is the orbital separation, and the 
orbits are assumed to be circular. For each BPS realiza- 
tion we have used Hurley's rapid binary evolution code 
(Hurley et al. 2002) to follow the evolution of 10 7 sample 
binary systems from star formation to the formation of 
the CO WD + He star systems based on three evolu- 
tionary scenarios (i.e., the He star, EAGB and TPAGB 
channels; for details see Wang et al. 2009b). We then 
apply the calculations described in ^ and assume that if 
the initial parameters of a CO WD + He star system are 
located inside the relevant contour of Figure [2l a double- 
detonation explosion occurs. 



log P i (days) 

Figure 2. Regions in the initial orbital period-secondary mass 
plane (logP 1 , Mr,) for WD + He star binaries that produce SNe la 
for various initial WD masses. 



is large, which leads to the rapid shrinking of the orbital 
separation. After about 11 million years, the helium star 
begins to fill its Roche lobe while it is still in the he- 
lium core-burning stage. The mass-transfer rate is sta- 
ble and with a low rate of ~2 x 10~ 8 Af©yr _1 , which 
results in the formation of a helium layer on the surface 
of the CO WD. After about 5 million years, the mass of 
the helium layer increases to 0.1 Af Q , at which point a 
detonation is assumed to occur at the base of the he- 
lium layer, and further assumed to produce a double- 
detonation explosion. At this moment, the mass of the 
helium star is A/f N = 0.50 Af Q and the orbital period is 
log(P SN /day) = -1.64. 

Figure 2 shows the initial contours for producing SNe 
la in the logP 1 — M\ plane for various WD masses, i.e., 
M WD = 0.8, 0.9, 1.0 and 1.1 Af , where P 1 and M\ are 
the initial orbital period and the initial mass of the he- 
lium donor star, respectively. From this figure, we can 
see that the contours are slightly shifted to higher periods 
with the increase of initial WD masses. This is because a 
helium star with a specific mass has a smaller Roche-lobe 
radius with a more massive WD companion. The left 
boundaries of the contours are set by the condition that 
RLOF starts when the secondary is on the helium ZAMS, 
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Figure 4. Delay time distributions of SNe la. The filled squares 
and open circles are taken from Totani et al. (2008) and Maoz et 
al. (2011), respectively. 

3.1. Rates and delay times 

We present results for a constant star-formation rate 
(SFR) - assumed to be constant for the past 14 Gyr - and 
also for a delta- function SFR, i.e., a single instantaneous 
starburst. In Figure 3, we show the evolution of Galactic 
SN la rates by adopting Z = 0.02 and SFR = 5 M^yr" 1 . 
This gives a Galactic event rate of ~0.6 — 1.8 x 10 yr" 1 , 
which is roughly one third of the inferred Galactic SN la 
rate (3 - 4 x 10 -3 yr" 1 ; Cappellaro & Turatto 1997). 
This is in good agreement with the measured rates of 
SNe lax (31±^ SNe lax for every 100 SNe la in a given 
volume; Foley et al. 2012). The birthrate from a cc A = 
0.5 is lower than that of a cc X = 1.5, since the post- 
CE binaries are more likely to be located in the SN la 
production region for a cc X = 1.5. As expected, these 
rates are consistent with the double-detonation model 
only producing part of the overall SN la rate (for a recent 
review of other potential SN la formation channels see 
Wang & Han 2012). 

The delay times of SNe la are defined as the time in- 
terval between the star formation and SN explosion. The 
various progenitor models can be examined by compar- 
ing the delay time distributions with that of observa- 
tions. Figure 4 displays the delay time distributions of 
SNe la for the double-detonation model. In the figure, 
we see that SN la explosions occur between «70 Myr and 




Figure 5. Distribution of peak bolometric magnitudes from the 
population model (for two different a ce A values). The data are ob- 
served SN lax V-band magnitudes from Foley et al. (2012), several 
of which are lower limits. 

«710Myr after the starburst, consistent with the current 
SNe lax sample, all of which have originated in late- type 
- i.e., star-forming - galaxies. 

3.2. Luminosity distribution 

To quantify the relationship between the WD explo- 
sion mass (A/wd) and the SN la peak bolometric mag- 
nitude (Mboi)j Ruiter et al. (2012) recently carried out 
a series of simulations for a range of WD masses based 
on ID sub-Chandrasekhar mass pure detonation mod- 
els, following the earlier similar calculations by Sim et 
al. (2010). We note that these simulations assumed a 
central detonation, with no specified triggering mecha- 
nism; we assume that their broad results are applicable 
to double-detonation events. The SN la peak brightness 
in their simulations is directly related to the WD explo- 
sion mass. We apply the A%d — M\, Q \ relationship found 
by Ruiter et al. (2012) to derive the SN la peak bright- 
ness. The resulting luminosity distribution is plotted in 
Figure [5] There we also show observed peak V-band 
magnitudes from Foley et al. (2012), although several of 
those are lower limits and the sample is still probably too 
small to draw strong conclusions from the shape of the 
distribution. Nonetheless, the predicted peak magnitude 
range reproduces the full observed diversity of SNe lax 
(-14.2 > M Y > -18.9 mag; Foley et al. 2012). The pre- 
dicted range (—13 > Mboi > —19 mag) does add a tail 
stretching to fainter magnitudes, which might indicate a 
deficiency in the models or an observational bias against 
discovering the fainter SNe lax. 

4. DISCUSSION AND CONCLUSIONS 

The most relevant known system to our current work 
is perhaps CD— 3011223, which has been identified as 
a CO WD + sdB star system with a ~1.2h orbital pe- 
riod (Geier et al. 2012; Vermes et al. 2012). Due to the 
short orbital period, angular momentum loss from GWR 
is large. After about 25 million years, the sdB star will 
begin to fill its Roche lobe while it is still in the He-core 
burning stage. CD— 3011223 may explode as an SN lax 
via the double-detonation model in the future evolution. 
Other known CO WD + He-donor systems (e.g., KPD 
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1930+2752, V445 Pup, and HD 49798 with its WD com- 
panion) may produce SNe la via other progenitor models 
(see Wang & Han 2012). 

The surviving companions from the model presented 
here for SNe lax may also explain hypcrvclocity helium- 
rich stars like US 708 (e.g., Hirsch et al. 2005) due to 
the short orbital periods at the moment of SN explosion 
(see also Justham et al. 2009; Wang et al. 2009). Study- 
ing such high-velocity helium-rich stars (and their WD 
descendants) might provide a way to test this model. 

In this Letter, we have systematically studied a he- 
lium double-detonation model for SNe lax in which the 
donor star is non-degenerate. (1) This model can natu- 
rally produce SN explosions with helium lines and with- 
out hydrogen lines. (2) The event rate agrees with the 
inferred rates of SNe lax. (3) The explosions in this 
model occur between «70Myr and «710Myr after the 
starburst, i.e., in relatively young stellar populations, 
consistent with the host galaxy morphologies of known 
SNe lax. (4) By adopting an existing relationship be- 
tween CO mass and SN luminosity for pure detona- 
tions, we find a SN luminosity range from —13 mag to 
— 19 mag, which compares well to the current diversity 
of SNe lax (-14.2 > My > -18.9 mag). The overall 
population properties resulting from this model therefore 
seem promisingly consistent with the known collection of 
SNe lax. 

Confirmation that helium doublc-dctonation explo- 
sions can produce a form of SN la would finally an- 
swer a fascinating astrophysical question, and help to 
constrain our understanding of explosion physics. How- 
ever, our results suggest a new problem: why would sys- 
tems with non-degenerate helium donors make SNe lax 
but not - or only rarely - those with He WD donors? 
Degenerate donors tend to produce higher mass-transfer 
rates than our non-degenerate donors (see, e.g., Han & 
Wcbbink 1999), and higher mass-transfer rates require 
less massive helium layers for ignition (see, e.g., Shen & 
Bildstcn 2009). So perhaps helium layer detonations in 
systems with degenerate donors typically involve helium 
layers which are not massive enough lead to a double- 
detonation? Alternatively, if double-detonations arc trig- 
gered in systems with lower-mass helium layers then per- 
haps the explosion properties are not sufficiently altered 
for them to be identified as peculiar? We intend to study 
this in future by adopting more sophisticated assump- 
tions about the properties of the helium layer at ignition. 
The 2002cx-like subclass of SNe la certainly deserves fur- 
ther detailed study. 
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